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1. Introduction 

In the last decade there have been several 
claims for high energy 7— ray emission between 
250 GeV and 0.75 PeV from the X-ray binary 
pulsar Hercules X-1, some pulsed at the X-ray 
period [1] [2] [3], others at a period significantly 
blue-shifted by an amount agreed upon by several 
observers [4] [5] [6]. Recently, there have been 
several detections at a different blue-shifted pe- 
riod by another group [7]. There are reports of 
long term DC emission at a somewhat higher en- 
ergy [8] [9], but these are not statistically very 
significant, and not corroborated by recent more 
sensitive searches [10] [11]. The possibility that 
this and other such sources (e.g. Cyg X-3) are 
variable on timescales of several years is not sur- 
prising, given the evidence of the irregularity of 
the intrinsic power source in Her X-1, on a sim- 
ilar timescale, deduced from optical observations 
spanning the last century [12] [13]. However, here 
we will be concerned with the convincing pulsed 
behaviour at the slightly lower energy, which is 
characterised by episodic bursts of emission last- 
ing typically 30 mins. 

The progenitors of the 7-rays are usually as- 
sumed to be protons, because electrons would ra- 
diate their energy much too quickly in the mag- 
netic field of the neutron star. In this case the 7- 
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rays are thought to be produced by interaction of 
the protons in some gas target in the system, fol- 
lowing the decay of any tt" produced in collisions. 
If this picture is correct, then we run into prob- 
lems with the huge implied proton luminosity, 
which appears to be higher than the total power 
available from accretion. [14]. For example, the 
fiux from the Cygnus Array [6] implies a 7-ray lu- 
minosity of 1.6 X lO^^erg s~^, if produced isotrop- 
ically. The parent proton luminosity is probably 
an order of magnitude higher [15], implying a pro- 
ton luminosity about 100 times higher than the 
estimated bolometric luminosity [16]. This en- 
ergy budget problem is best balanced by invoking 
a narrow beam, which also has the advantage of 
explaining the episodic nature of the emission [14] 
if it wanders around in direction, possibly driven 
by the precession of the neutron star. 

There are also problems in explaining the high 
energies that individual protons appear to be 
achieving. Observation of 7-rays up to energies of 
750 TeV implies the acceleration of protons to an 
energy several times higher [15], probably several 
PeV. There have been several models proposed to 
explain the acceleration of protons to energies of a 
few PeV, but in general they have been shown to 
be untenable in the uniquely well understood en- 
vironment of Her X-1. Shock acceleration above 
the pole of the neutron star [17] [18] fails because 
of the enormous magnetic field which inhibits the 
formation of scattering centres [19] [20]; the pul- 
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sar wind shock mechanism [21] does not work be- 
cause the pulsar spins too slowly [20]; the disk 
dynamo [22] [23] becomes very inelEcient and the 
maximum energy available is much reduced for 
the high magnetic field found here [20]; turbulent 
reconnection of magnetic field may accelerate par- 
ticles in X-ray binaries [24], but the total amount 
of energy available for particle acceleration is ex- 
tremely low, and cannot explain the duration of 
the observed bursts [25]: the total amount of en- 
ergy associated with the magnetic field ouside the 
Alfven radius (where the plasma can still distort 
and wind-up the field) is much less than the en- 
ergy requirement inferred from the high energy 
7— ray observations. 

A pulsar dynamo model [26] successful in ex- 
plaining 7-ray emission from rapidly spinning iso- 
lated pulsars has been extended to slower, accret- 
ing pulsars [27], but there may be problems in 
maintaining vacuum gaps in the presence of ac- 
creting plasma and the hot plasma of the accre- 
tion disk corona which is likely to form [28] [29]. 

The interaction between cosmic-ray protons 
and Alfven waves has long been known to play 
an important role in cosmic-ray transport in the 
galaxy. Under certain circumstances, cosmic-rays 
generate Alfven waves, and their velocity is lim- 
ited to that of the waves. In other situations 
the cosmic-rays can be efficiently scattered by 
Alfven waves, and in the presence of both forward 
and backward waves, Fermi acceleration can oc- 
cur [30]. Such acceleration by fast-moving Alfven 
waves in the magnetosphere of accreting neutron 
stars is a possible method of achieving proton 
energies of 1 PeV [31], and here we investigate 
the application of this mechanism to the accre- 
tion column of Her X-1. The preliminary results 
of this work have been described earlier in less 
detail [32]. The fiuctuations which propagate as 
waves can be caused by the accretion splash on 
the neutron star, and travel outwards. Instabil- 
ities at the Alfven radius can also create distur- 
bances which can travel inwards. Refiection of 
protons by the converging field lines provides an 
opportunity for additional energy gain, as mir- 
roring traps the particle in the system for longer, 
together with a focussing mechanism by which 
the protons emerge in a narrow range of angles. 



The advantages of acceleration in the accretion 
column itself have been discussed elsewhere [14]. 
Since this work was first reported. Smith, Katz 
and Diamond [33] have extended their work to 
accretion columns in general, and although the 
approach is different from that given here, the 
results regarding maximum energy are in agree- 
ment. 

2. The Environment 

Her X-1 is an accreting X-ray binary pulsar 
whose environment is well understood. It lies 5.8 
kpc distant [12], and consists of a neutron star 
spinning with a period P„ = 1.24 s, which orbits 
its companion every 1.7 d [34]. The neutron star 
is surrounded by a tilted and precessing accretion 
disk [35]. Her X-1 is unusual in that the mass of 
the primary is well measured, M„ = 1.3M0 [36], 
and is consistent with that expected for neutron 
stars [37] [38], for which a radius of R„ = 10 km 
is typical. Cyclotron lines in the X-ray spec- 
trum indicate the presence of a magnetic field 
at the pole of the neutron star with fiux density 
_Bq = 5 X 10^^ G [39], in line with expectation 
for such an object. Inside the Alfven radius, Ra, 
which is coincident with the co-rotation radius at 
~ 200 Rn, this strong magnetic field controls the 
fiow of gas. The gas slides as an ionised spray 
[40], its geometry described by the magnetic field 
lines, and forms an accretion column, landing on 
the surface of the neutron star over an area of 
somewhat less than 1 km^. In this investigation a 
simplified version of the accretion column is used 
as follows: Let r be the distance from the cen- 
tre of the neutron star, and x the distance from 
the column centre. The column extends between 
r = R„ to r = Ra with circular cross sectional 
area A = Ao{r/Rn)^, radius Xc = (A/irY/^, free 
fall gas velocity Vg = Vo{r/Rn)~^^^ and magnetic 
field strength B = Bo(r/ R„)~^ . The zero sub- 
script refers to the value at the neutron star pole: 
5o = 5 X 10^2 Gauss, and Vo = ^/2GM„/R„ = 
1.3 X 10^*^ cm s~^ are known. The density at the 
pole is estimated from equipartition between the 
magnetic pressure and the ram pressure of the 
accreting gas at the Alfven radius, and then ap- 
plying continuity to give po = O.Olg cm~^. The 
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area of the column at the pole can be estimated 
from the disk thickness at the Alfven radius to be 
Ao = 0.03 km^ [25]. These figures yield a mass 
fiux M = 4 X 10^'^ g s~-^, and a total available 
gravitational power of 7 x 10^® erg s""'^ per pole, 
which is consistent with the bolometric luminos- 
ity, 2 X 10^^ erg s""'^ [16]. The temperature of gas 
required to produce the observed luminosity from 
the polar area is 1.5 x lO^K, in agreement with 
that measured [34]. 

Since the gas in the accretion column is ex- 
pected to be fully ionised, Alfven waves will prop- 
agate with a velocity va given by 



(1) 



after [41]. The velocity of the Alfven waves is 
usually much greater than the velocity of the in- 
falling gas, and so the latter is treated as station- 
ary in this investigation. The quantities p, B and 
Va are assumed independent of the distance from 
the column centre, x. 

2.1. Spectrum of Alfven Waves 

Two modes of Alfven waves are considered: one 
mode moving outward in the -\-r direction and one 
moving inward in the — r direction; these will be 
referred to as the '-|-' and '-' modes, and both 
have the same va = VAir). The amount of power 
passing through the column in the form of the 
two modes of Alfven waves is denoted and 
Here we assume that the '-|-' mode origi- 
nates at the neutron star surface, the waves in 
the field perhaps caused by the accretion splash. 
The smallest wavenumber for this mode is taken 
to be ko+ = 10/ Rn- The '-' mode is assumed 
to start at the Alfven radius, the waves possi- 
bly initiated by the Rayleigh- Taylor instabilities 
at the interface between orbital motion in the 
disk and free fall motion in the accretion column 
[42]. The smallest wavenumber for this mode is 
taken to be kg- = 10/i?^. In other words, the 
longest wavelength in each spectrum is approxi- 
mately one tenth the size of the region causing 
the disturbances in the field lines. These minu- 
mum wavenumbers are uncertain, and are varied 
to assess the effect on the acceleration process. 

To calculate how strongly a particle of a given 
energy is scattered, we need to know how the en- 



ergy density is distributed over the spectrum of 
Alfven waves. We will assume that the spectrum 
takes the form 
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where h is the amplitude of the waves at a given 
wavenumber, k. Hence we can find the total en- 
ergy density resident in the spectrum, U, by in- 
tegrating between ki and k2 which, if k2 » ki 
and n < —I reduces to 



U 



n + l 
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Now if we specify that the total magnetic energy 
fiux in Alfven waves to be a constant L at any 
distance from the neutron star, then since L = 
UvaA we have 
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for any assumed L, n and ki. We also have to 
allow for the fact that the velocity of the wave, 
Va changes with r, and hence a wave with initial 
ki = kg and va = va^o (at the source of the spec- 
trum) will become ki = koVA^o/vA as va changes. 
Finally we have 
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which from equation 2 gives us d(h'^)/dk for each 
mode as a function of va, once L, n and kg are 
specified, A having being found in section 2. Of 
the variables, the total power in the Alfven waves, 
L, is the hardest to determine. 

2.2. Electric and Magnetic Fields Associ- 
ated with Alfven waves 

Consider an Alfven wave propagating in a 
plasma along the direction of a strong ambient 
field. Seen in the frame in which the plasma 
has no bulk motion, the total magnetic field is 
(B -|- b), where b is an oscillating component or- 
thogonal to the ambient field B. If the conductiv- 
ity is very high, there is also an oscillating electric 
field E which is in phase with b, and orthogonal 
to B and b, travelling with a velocity va such 
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Figure 1. The Lorentz factor of the Alfven waves, 
7^, as a function of the distance from the neutron 
star, r, in the accretion column. The distance r 
is given in units of the neutron star radius, 



that E = VAh- The velocity of the wave rela- 
tive to the plasma is given by equation 1 (At low 
speeds, when p/B'^ is high enough, this expres- 
sion reduces to the familiar B /(iirpy/^). The 
electric and magnetic field may be found for any 
other frame moving at a velocity v along B by 
using the Lorentz transforms. In the rest frame 
of the wave, moving at a velocity va along B, the 
transverse fields are 

E' = jA(E-VAb) = (6) 

and 

b' = jAib-VAE/c''), (7) 

where is ja is the Lorentz factor of the waves. 
Similarly the inverse transforms become 

E = jAVAb' (8) 

b = 7Ab'. (9) 

The behaviour of ja in the accretion column 
is shown in Fig. 1. It can be seen that in the 
inner half of the accretion column, the waves be- 
come relativistic. In the extreme limit, the elec- 
tric field associated with the wave is E = cb, and 



the Alfven wave takes on the form of an electro- 
magnetic wave in free space. The circumstances 
under which this transition has occurred are un- 
usual: the condition B^ / p >> AttP, obtained 
from Eqn. 1, is normally satisfied because p ^ 0, 
but here the density is actually increasing. It is 
the rapid rise of B which causes the transition in 
this case. 

3. Calculation of Scattering 

The scattering of protons by Alfven waves de- 
pends on the total force on the particle. For 
slowly moving Alfven waves, the magnetic field 
only need be considered. For relativistic waves, 
the electric field becomes appreciable, and the 
magnetic and electric forces become comparable. 
In this investigation the scattering is calculated 
for the particle in the frame of the wave. Here, 
the electric field E' is zero, and only the mag- 
netic field b' remains. Once the magnetic field 
in the frame of the column, b, is known, that in 
the frame of the wave, &', is found from equation 
9. The calculation of the scattering in this frame 
then proceeds as described in the next section. 

3.1. Pitch Angle Scattering Due to Alfven 
Waves 

Consider a particle spiralling in a stationary 
magnetic field of strength B, so that the parti- 
cle has a Larmor radius = p/eB, pitch angle 
6, and hence pitch equal to 27rriCos6'. Suppose 
a spectrum of small scale tranverse waves is su- 
perposed. At a particular time, at the location of 
the particle, this transverse field will be of magni- 
tude b. In a time At the particle moves a distance 
s = cAt, and the pitch angle changes by 

Ae = —-. (10) 

vlB ^ ^ 

The total effect on the particle depends strongly 
on wavelength. A: for relatively short wave- 
lengths, such that A << r^, the waves add in- 
coherently to give a very erratic transverse field. 
We can assume defiections are uncorrelated over 
even quite short distances, and the effective field, 
bgff oc ^/s. At the other extreme of wavelength, 
when A >> r^, b does not change appreciably as 
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the particle turns round its helix. Effectively the 
total field is changed gradually, without an ap- 
preciable change in magnitude, and the guiding 
centre will follow the disturbed field line adiabat- 
ically. Here we have no change in pitch angle, very 
little displacement (unless very long wavlength 
waves cause field line wandering) and hence no 
cross-field diffusion. 

However, when A ~ r^, the particle sees a con- 
stant transverse field all the way round its orbit, 
which causes a large change in pitch angle. If we 
decompose the wave into right and left circularly 
polarised components, one set will have a radial 
component of h which turns with the particle and 
so gives a constant d9/dt, and hence A6' cx s in- 
stead of ~/s. In fact, integrating over the whole 
spectrum 



Mr 



B 



(11) 



such that A6' cx for scattering by waves where 
A < and zero for longer wavelengths. The 
effective field, h^jj is given by 



V dk 



(12) 



k=K 



where K = l/27rriCos6' is the inverse of the 
pitch, or the effective wavenumber of the parti- 
cle. These last two equations are equivalent to 
the treatment of [30] which was more hydrody- 
namical in nature. 

The cross-field displacement is also important 
when A < r^, and in fact the guiding centre is 
displaced in a random azimuthal direction by an 
amount 



dB _ ^dr 
~B ~ 7' 



(14) 



From the adiabatic invariant sin^6' cx B, if the 
proton travels a distance s we have by differenta- 
tion with respect to s 



1 d(sm^0) sin^e dB 







B ds 52 ds 

whence from dr/ds = — cos 6 and equation 14 

d9 3 sin 9 
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ds 2 r 
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Of course, no displacement of the guiding centre 
occurs in this case. 

3.3. Transformations and particle energy 
gain 

The mechanism for energy gain is similar 
to second-order Fermi acceleration, except that 
here the scattering centres are relativistic Alfven 
waves, which provides an opportunity for very 
rapid acceleration. An Alfven wave of amplitude 
h travelling with velocity va and Lorentz factor 
jA in the accretion column frame has amplitude 
h' = h/jA in the frame of the Alfven waves. In 
transforming to the wave frame a particle with 
energy velocity /3c and initial longitudinal mo- 
mentum Pi in the column frame becomes 



Pi = lAipi - I3a£z) 



(17) 
(18) 



in the frame of the wave (quantities primed). Af- 
ter accounting for the scatter of the particle we 



r^AO sin ( 



(13) 



This is necessary to determine whether a particle 
wanders sideways out of the accretion column. 

3.2. Pitch Angle Change Due to a Con- 
verging Field 

The pitch angle is also changed, this time sys- 
tematically, by the converging magnetic field of 
the neutron star. If the proton is near the mag- 
netic axis at a distance r from the neutron star. 



then B 



and hence 



PJ = IaIPj +PA£'i) 
£f = 7a(£- +l3AP'f) 



(19) 
(20) 



to transform back to the column frame, using the 
subscripts i and / to denote the initial and fi- 
nal states. The energy remains unchanged in the 
frame of the wave, e[ = e'j, and the change in p' 
depends on the change in pitch angle. In fact 
for total refiection, where p'j = —p[, Ae/e = 
2'^\(3a{P -\- Pa), which reduces to Ae/e = 47^ 
for relativistic protons and waves. From the high 
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values of ja attained in the accretion column as 
indicated in Fig. 1, it is evident that we might 
expect large energy gains: ja ~ 100 implies 
As/s ~ 4 X lO'*. Such energy gains can only be^T" 
achieved if the irregularities in the field are strong > 
enough to be able to turn the particle round in 
the space available, and if the particle is not to ^ 
wander sideways outside of the column. ^ 

■a 

4. The Simulation 

Protons are injected isotropically into the ac- 
cretion column of section 2 with initial energy 
£j and position Tj-. A transformation is then 
made into the frame of one of the two Alfven 
wave modes. Here the particle is moved by a 
small amount s, and the pitch angle scattering 
and sideways displacement are calculated. The 
amount of pitch angle scattering is determined by 
the method of rejection, by sampling a gaussian 
distribution with width given by 6rms in equation 
11. This displacement is optimised to correspond 
to a reasonable pitch angle change so that the 
particle is not followed to undue precision. The 
contribution to the scattering from the other wave 
mode is then calculated, and then finally, back in 
the frame of the plasma, the adiabatic pitch angle 
change is accounted for. This process is repeated 
whilst keeping track of the particle position and 
momentum, until the particle either leaks from 
the side of the column (x > Xg), hits the star 
(r < R„) or escapes (r > Ra)- Many protons 
are thus injected to build up a spectrum of final 
energies for each of the three possible fates. 

5. Results 

5.1. The Canonical Spectrum 

The canonical parameters for this investigation 
are those mentioned in section 2 together with 
n = -3/2, L+ = L_ = 10^5 ergs"!, r,- = Ra/I 
and £j =1.1 GeV. The power in the two modes 
of Alfven waves is set at roughly 1% of the to- 
tal gravitational power available per pole (_Lj/2), 
which at the Alfven radius, Ra, corresponds to an 
energy density U c:^ Ub, the background magnetic 
energy density, and hence represents a realistic 
upper limit to the power available. At the pole 
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Figure 2. The spectrum of protons which emerge 
at the Alfven radius after acceleration in the ac- 
cretion column. 



of course, the energy density of the waves is far 
less than that in the ambient field, U ~ 10~^Ub- 
Concerning the other parameters, the spectral in- 
dex of the magnetic irregularities corresponds to 
a Kolmogorov spectrum and the initial energy of 
the particles is approximately equal to the energy 
that the protons achieve in free-fall. 

The proton spectrum per injected proton, 
dn/de, obtained for these parameters is shown in 
Fig. 2. The majority of protons in fact emerge 
from the column at Ra , only a small fraction leak- 
ing sideways and very few hitting the star. Most 
protons emerge with more than 10 TeV of en- 
ergy, and the mean energy obtained is 30 TeV. 
The spectrum extends up to around 250 TeV with 
a slope of -2.0, after which point it steepens to 
around -3.0. Above this steepening there is no 
evidence for an absolute cut off. The statistics 
limit the maximum proton energy attained in this 
investigation to just under 1 PeV. 

The energy and position of a typical proton are 
shown as a function of time in Fig. 3. This par- 
ticular proton achieves an energy of 30 TeV in 
10 ms. 
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Figure 3. Some features of the history of a typical 
particle are shown. Both (a) the energy, e, and 
(b) the distance from the neutron star r are given 
as a function of time since injection, t. This parti- 
cle escapes after 10 ms with an energy of 30 TeV. 
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Figure 4. The mean escape energy, e, versus the 
particle injection distance from the neutron star. 
The distance is given in units of the neutron star 
radius, The dashed line indicates the canon- 
ical value used in the investigation. 



5.2. Sensitivity of spectrum to input pa- 
rameters 

The parameters Bg, Ag, Po and Vg are assumed 
to be sufficiently well known that they are not 
altered in this investigation. The remainder are 
not well known, however, and each is varied in 
turn whilst holding the others at the canonical 
values. The emergent spectrum is obtained for a 
range of Tj-, n, ko+ and kg-, -t+ and 

The injection point, ri, was varied between 
2 R„, near to the neutron star, out to near the 
Alfven radius, Ra, and the corresponding mean 
escape energy is given in Fig. 4. The mean escape 
energy shows some dependence on Tj-: the energy 
is degraded markedly as ri approaches Ra, but 
this is to be expected as the protons tend to es- 
cape before they have been significantly acceler- 
ated. The mean particle energy is enhanced to 
200 TeV as the injection point moves a closer to 
the star than the canonical value, i?A/2. There- 
fore, if the protons had been injected uniformly 
along the column, the mean energy would have 
been nearer this value. 

A range of n from -1.0625 to -2 produced no 
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Figure 5. The mean escape energy, e, versus the 
minimum wavenumber at the source of the Alfven 
waves, A;o_|_(full) and kg- (dotted), both given is 
a fraction of the canonical values described in the 
text. 



Figure 6. A contour map of the mean escape en- 
ergy, e, as a function of the power in both the in 
(-) and out (+) wave modes. The power is ex- 
pressed as a fraction of the accretion luminosity 
per pole (_Lj/2). The contour lines are labelled in 
units of TeV. 



dramatic changes in the spectrum. The spectrum 
is slightly degraded for the most negative values 
of n. More important seems to be the power con- 
tained in the Alfven waves, and rather 
than the shape of the spectrum. The range of 
power tried for these parameters spans 10^® to 
10^'^W in view of the uncertainty in the mech- 
anism which mechanically couples the inflowing 
gas to the field lines to produce Alfven waves. 
The mean energy obtained as a function of the 
power residing in the two Alfven waves modes is 
illustrated in Fig. 6. The spectrum is not en- 
hanced significantly when both and L- are 
above the canonical values. However, as L- drops 
there is a degradation in energy, and a tendency 
for more protons to leak from the column before 
emerging at Ra- If i+ is lowered, the mean en- 
ergy of the protons drops somewhat faster. How- 
ever, the energy degradation is severe only for 
very small values of L: remarkably, even with 
just lO^s W in each of the modes (10-*i.j/2), the 
protons still emerge with a mean energy of nearly 
3 TeV. 

The cut-off wavelength present in the spectrum 



of Alfven waves is hard to estimate, and so ko+ 
and kg- are tried at factors of up to 100 either 
side of the canonical values, and the results given 
in Fig. 5. The spectrum is insensitive to kg-, 
but the maximum achieved energy goes up dra- 
matically as ko+ is lowered. This effect might be 
expected: the maximum attainable energy for a 
particle is achieved when the gyroradius of the 
particle is comparable to the longest wavelength 
present in the spectrum of the Alfven wave mode 
by which it is being scattered. For the '-|-' mode, 
the canonical value ko+ is equivalent to a wave- 
length of 0.1i?„, which is the same as the gyrora- 
dius of a particle of energy 40 TeV at i?A/2. This 
is just below the break in the spectrum of acceler- 
ated protons. Particles which undergo scattering 
a little further in towards the star can reach the 
higher energies seen in the simulation. If ko+ is 
decreased then this maximum energy is seen to 
rise as expected. For maximum wavelengths of 
1 Rn in the '-|-' mode, the mean energy is en- 
hanced to above 100 TeV. In contrast, the wave- 
length associated with the incoming waves is al- 
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Figure 7. The distribution over solid angle, 
dn/dil for those protons which emerge at the 
Alfven radius. The data is expressed as a frac- 
tion of that expected for an isotropic distribution 
(dashed line). 



ways large enough to interact with protons in the 
energy range of interest, and the mean energy 
of accelerated particles shows almost no change 
when this parameter is varied by a factor of a 100 
either side of the canonical value. 

The funnel-shape of the magnetic field in the 
accretion column has the important effect of sys- 
tematically reducing the pitch angle of the par- 
ticle through adiabatic invariance, which keeps it 
in the system for longer by inhibiting leakage out 
of the side of the column. More importantly, by 
this same process the protons are also focussed 
into a beam: by the time they emerge at Ra, 
the protons are travelling roughly parallel to the 
magnetic field, as shown in Fig. 7, where the dis- 
tribution in solid angle is given, dN/dil together 
with that expected for an isotropic distribution. 
For those particles travelling with zero pitch angle 
along the accretion column the intensity is over 
17 times higher than if the protons had emerged 
isotropically. The particle beam would stand the 
highest chance of being detected when directed 
at Earth, presumably following the production of 
7— rays after passing through gas in the accre- 



Figure 8. The distribution of acceleration times, 
tacc, for those protons which emerge at the Alfven 
radius. 



tion disk. In this case, detection would therefore 
cause an overestimation of the power of the ac- 
celerator by this substantial factor, if the usual 
assumptions of isotropic emission are made. 

For the canonical spectrum, the time spent in 
the system for each proton is shown in Fig. 8. 
Protons leave the system in their accelerated state 
after typically 10 ms. The process is therefore 
very fast compared to the timescale of the neutron 
star rotation period 1.24s, and significantly faster 
then the infall time of the accreting plasma from 
the Alfven radius to the neutron star, 

[got:) ^^^^ 

which is approximately 90 ms. 

6. Loss Mechanisms 

In the simulation, no energy losses were as- 
sumed. Here we check various possible meth- 
ods of energy loss which could restrict the max- 
imum energy attained. The protons are subject 
to energy losses by synchrotron radiation in the 
strong magnetic field of the neutron star, interac- 
tion with X-rays from the hot-spot, and inelastic 
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collisions with the gas in the accretion column. 

6.1. Losses through pp collisions in accre- 
tion column 

The density of gas in the accretion column, p is 
given in section 2. A proton travelling from a dis- 
tance r from the neutron star, outwards through 
the accretion column will traverse an amount of 
matter S given by 

/ \ -3/2 
E = 6600 (-^) g cm-2 (22) 

Protons travelling from a point r/R„ = 20 axi- 
ally outwards will thus traverse 75 gcm~^ of mat- 
ter. However, the protons are accelerated in a 
region further out than this, and they may en- 
counter substantially less than this amount. On 
the other hand, the protons do not travel directly 
out of the system. In fact, from the simulation, 
the distribution of the amount of matter traversed 
is shown in Fig. 9. Typically less than 15 g cm~^ 
is encountered (which indicates that the protons 
get no closer than about 60 i?„ to the neutron 
star). This is considerably less than the pp inter- 
action mean free path, ~ 80 g cm~^ and so the 
accelerator will not suffer serious losses through 
this channel. However, this amount of matter 
may serve as a thin target for continuous weak 
7-ray production. 

6.2. Synchrotron Loss 

The synchrotron loss time is given by 

= 2^5^ = ^ lO-'^it^s"' s (23) 

where the proton is travelling orthogonal to the 
magnetic field [43]. If we define the mean free 
path for synchrotron loss as = ctg , expressing 
it as a fraction of the distance from the neutron 
star, we get 

assuming a dipole field, B oc r~^, and an average 
pitch angle of 7r/4. An indication of serious en- 
ergy loss is given when A^/r is less than unity i.e. 
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Figure 9. The distribution of the total grammage, 
S, encountered by those protons which emerge at 
the Alfven radius. 



the mean free path is less than the scale size of the 
region. This occurs for r < 12R„ for 1 TeV pro- 
tons, and r < 50i?„ for 1 PeV protons. Since the 
acceleration generally occurs in a region further 
out than this, and most of the time the protons 
are travelling with very small pitch angles, energy 
loss by this process is not important except pos- 
sibly for the very highest energy protons when 
they get relatively close to the neutron star, in 
which case we might expect secondary 7-rays to 
be produced. In this way synchrotron loss could 
ultimately limit the maximum energy achievable 
by this acceleration mechanism to around 1 PeV. 

6.3. Losses through photopion production 

Photopion production has been shown [44] to 
be a possible source of 7-rays for protons above 
1 PeV in X-ray binaries. The reaction p + j ^ 
p + TTo has a threshold given by 

£p£x > 1-5 X 10^^(1 - cos a)-^ eV^ (25) 

where Sp and ex are the energy of the proton and 
photon respectively in eV, and a is the angle be- 
tween them. Since the temperature of the hotspot 
on the neutron star is 10® K, the average energy 
of the emitted photons is < ex >— 25 keV, and 
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so the threshold energy for this reaction is then 
3 TeV for head on collisions, and rises to 400 TeV 
when a ~ 10°. As the protons are bouncing back 
and forth in the accretion column they will meet 
the X-ray flux at the full range of angles, and so 
this reaction is in the energy region of interest. 
Above this threshold the mean free path of a pro- 
ton is given approximately by 



An 



600 



L 



X 



]^Q37gj.g g-1 



ex cm (26) 



where is the distance from the hot-spot (rather 
than the distance from the centre of the neutron 
star). Substituting for ex and the X-ray luminos- 
ity Lx ~ -Lj = 2 X 10^^ erg s"-'^, and expressing 
as a fraction of the radius at that point. 



r ■ \Rr. 



(27) 



which is greater than unity except very close to 
the neutron star. 

In addition, the protons will be screened from 
the X-ray flux by the gas in the column. The 
semi-thickness of the column, S^, may be ex- 
pressed as 



p{r)xc{r) 



100 1 — 



gem 



(28) 



which is greater then the Thompson value of 
2.5 gcm~^ for the inner part of the column. As- 
suming flrstly that the protons are concentrated 
in the centre of the accretion column, and sec- 
ondly that the X-rays impingeat an angle normal 
to the column, the overall interaction length re- 
duces to 



7.5r exp 



40 



Rn 



(29) 



which is signiflcantly greater than unity for the 
whole length of the column. This is illustrated in 
Fig. 10, together with the interaction length in 
the case of no screening. The protons will prob- 
ably be distributed more evenly in the column of 
course, but this will make little difference at large 
distances since even with no screening the inter- 
action length there is still very large. Even if the 
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Figure 10. The p7 interaction length as a 
function of the distance from the neutron star, r. 
The interaction length is expressed as a fraction 
of r, and is shown with and without the screening 
effect of the gas in the accretion column (full and 
dashed lines respectively). 



column is curved overall, closer in it stands verti- 
cally above the hot spot, and the screening will be 
very much more effective because the X-rays will 
not be coming in the sides, but up through the 
lower, densest part of the column. For these rea- 
sons, the interaction length given in the Figure is 
probably conservative. Therefore, partly because 
of the screening effect of having the accelerator in 
the gas of the accretion column, signiflcant loss by 
photopion production does not occur. As a con- 
sequence of this, secondary 7-ray production by 
this mechanism will also be suppressed. Photo- 
pair production will also occur with a similar en- 
ergy loss rate, but at substantially lower energies 
than the pion production and so this channel is 
also unimportant. 



7. Estimate of 7-ray flux 

The flux of 7-rays from Her X-1 detected by 
the various experiments can be summarised as fol- 
lows. Around 1 TeV, long term fluxes of several 
10~"'^"'^cm~^s~"'^ and bursts up to 10~®cm~^s~"'^ 
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are seen, whilst at several hundred TeV, long term 
fluxes around 10~^^cm~^s~^ and bursts up to 
10~"'^"'^cm~^s~"'^ have been reported. We now con- 
sider the acceleration process described here in 
view of these fluxes. 

Signiflcant 7-ray production depends on the ex- 
istence of a sufficiently dense gas target to act as 
a beam dump, and the exact flux depends on the 
details of the hadronic cascade which will be initi- 
ated when the protons strike such a target. Here 
we will simply assume that a fraction ~ 10% 
of the proton luminosity is converted to 7-rays. 
The inner edge of the accretion disk may func- 
tion as the target, and most of the time may be 
too dense or too thin to act so efficiently, but cer- 
tainly some of the time this figure will serve as a 
good approximation. 

After passing through such a target, the flux of 
7-rays measured at Earth will be 

P^^(^^_jM^^j^~^E,) cm^-s-i (30) 

where Q is the enhancement due to proton beam- 
ing; tA and fp are the fraction of the bolometric 
luminosity converted to Alfven waves and then 
high energy protons repectively; E'p is the mean 
energy of high energy protons produced by the 
accelerator; /(> E'p) is the fraction of protons 
emerging with energy greater than E'p; and fl- 
nally, D is the distance to the source. Using 
Q = 17, = 10%, D = 5.8 kpc and Ep ~ 30 TeV 
we obtain 

= 1.7 X lO-^°£p£Af(> Ep) cm^-s-^ (31) 

Assuming TeV 7-rays are produced by 10 TeV 
protons, then since /(> lOTeV) ~ 1 we require 
^p^A > 0.1 to reproduce the long term low flux 
at this energy. Even with £p£A = 1, it is not 
possible to account for the burst flux. At the 
higher energy, 7-rays of 100 TeV require parent 
protons near 1000 TeV, where /(> lOOOTeT/) is 
only 10~^, and even the low fluxes reported are 
only explained if fpf^ = 1. 

Clearly, the observed luminosity of the 7-rays 
is still a problem. Either the beaming factor of 



the protons needs to be larger than found in this 
simulation, or the power available from accretion 
must occasionally be much greater than the aver- 
age bolometric luminosity. Since the latter is al- 
ready Eddingtonian, we require occasional super- 
Eddington accretion rates to explain the observa- 
tions. Such a situation is not so unrealistic since 
the accretion geometry is decidely asymmetric. 
(This requirement is true of course for any model 
which ultimately draws its power from the gravi- 
tational infall of matter, and the extent to which 
we must appeal to super-Eddingtonian accretion 
is much less in a model like this where the pro- 
tons are naturally beamed.) An implication of in- 
creased accretion rates is that the gas in the accre- 
tion column will be accordingly denser, which will 
improve the production of 7-rays in the column 
itself up to a point, but ultimately pp collisions 
will degrade the maximum energy achievable. 

8. Conclusions 

A simple model for the acceleration of pro- 
tons by scattering off relativistic Alfven waves has 
been set up in the well understood environment 
of the accretion column of Her X-1. It may be ap- 
plicable to other similar accreting X-ray pulsars 
and possibly some cataclysmic variables for which 
TeV emission has also been detected, but Her X- 
1 remains the prototype of this type of source 
and the only one for whose parameters are very 
well established. The mechanism is reminiscent 
of classical second-order Fermi acceleration, ex- 
cept that here the "magnetic clouds" are moving 
relativistically. When plausible parameters are 
assumed for the spectrum of Alfven waves and 
the initial state of the particles, most protons are 
accelerated to energies in excess of 10 TeV, with 
a mean energy of 30 TeV, and an spectrum 
which extends beyond 100 TeV before steepen- 
ing. Particle injection slightly nearer to the neu- 
tron star shifts the mean proton energy to around 
200 TeV, and uniform injection of protons along 
the column would raise the mean energy to this 
value. The maximum energy seems to be mainly 
limited by the maximum wavelength of the Alfven 
waves. The presence of waves with wavelength 
~ R„ enhances the mean energy to 250 TeV. 
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The highest individual energy achieved is over 1 
PeV, limited only by the statistical nature of the 
simulation, although the spectrum has steepened 
markedly by this point. Whilst mean energies of 
over 100 TeV are achieved in favourable circum- 
stances, 30 TeV seems achievable following rea- 
sonable assumptions about the environment, and 
production of 1 TeV protons seems to be an al- 
most inevitable by-product of the accretion mech- 
anism in this system. 

For the simulation, the single test-particle ap- 
proach has been followed. No attempt has been 
made, however, to consider the energetics of a 
flow of such protons: how they interact back on 
the system, and exactly what luminosity can be 
supported in such a flow are very complex ques- 
tions that require further work. Nor has the ques- 
tion of particle injection been addressed. In the 
simulation, relativistic protons were injected into 
the mid-point of the accretion column, but of 
course the injection could occur at any point or 
even uniformly along the column. The mecha- 
nism for particle injection is yet another question. 
Here we simply assume that some of the rela- 
tivistic particles produced in the accretion pro- 
cess flnd their way into the upper reaches of the 
column. 

The process has the advantage of producing a 
focussed beam of protons which substantially re- 
duces the power requirement of high energy par- 
ticles. The beam emerges in the direction of the 
inner accretion disk, where there is sufficient ma- 
terial for the production of 7-rays, and which may 
help explain the blue-shifted nature of the pulsed 
emission [45] . Explaining the flux of observed 7- 
rays is still a problem however, even with the 
beaming. A pulsed signal of other secondaries 
is possible, including neutrinos. Sporadic detec- 
tion of secondaries is easily explained if the beam 
wanders around for some reason e.g. due to the 
precession of the nuetron star and hence the mag- 
netic held. Constant low-level production of such 
secondaries is also possible, due to the thin gas 
embedded in the accretion column itself. 
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